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Abstract: we show that the solution to an oscillatory-singularly perturbed ordinary differential equation 
may be asymptotically expanded into a sum of oscillating terms. Each of those terms writes as an oscillating 
operator acting on the solution to a non oscillating ordinary differential equation with an oscillating correction 
added to it. 

The expression of the non oscillating ordinary differential equations are defined by a recurrence relation. 
We then apply this result to problems where charged particles are submitted to large magnetic field. 



1 Introduction and results 

Purpose 

The goal of this paper is to manage and justify the asymptotic two scale expansion as e — > of the 
solution X e (t) = X £ (i; x, s) to the following singularly perturbed dynamical system: 

^£ = a(t,— ,X £ ) + ib(i,X E ), X £ ( S ;x, S )=x, (1.1) 
at e e 

and to apply this to models of charge particles submitted to a strong magnetic field which is possibly 
non uniform. 



More precisely, C™ standing for the space of functions that have continuous and bounded deriva- 
tives until the order m, we assume that 

a( v ,0 € (C*™ +1 (M x R x R d )) d , 9 ^ a(t,6»,x) is 27r-periodic for every t € K and x E R d , (1.2) 

b( v )e(Cj m+2 (lxKY, (1-3) 

for some m > 0. In particular, it implies that \V X ■ b| < C, uniformly on 1 x if We also suppose 
that the solution Z(t, 6;z) to 

— =b(t,Z), Z(*,0;z)=z, (1.4) 

is known and 2ir— periodic in 9 for every t £ M and z £ M. d . 

Under those assumptions, we prove that X e (-;x, s) : R — > M. d admits the following expansion: 

X e (t; x, s) = X°(t, ^— s x, s) + eX 1 ^, t —^-; x, s) + e 2 X 2 (i, x, s) + . . . , (1.5) 
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as e — > 0, where the X l (t, 0; x, s) are 2n— periodic in 9. The initial condition in (jl.ip is assigned to 
X° by: 

X°(s,0;x,s) = x and X l (s,0;x,s) = for i > 1. (1.6) 

Motivations 

The target we have in mind while writing this paper is the conception of pic-methods for plas- 
mas submitted to large magnetic field. Direct applications of this are the simulation of magnetic 
confinement fusion or isotope resonant separation experiments. 

A plasma submitted to a strong magnetic field could be modelled by several Vlasov equations, 
according to experimental conditions. For instance if we are interested in describing the global 
behaviour of a plasma in a magnetic confinement fusion experiment, the following equation 

^ + v • V,f + ((E(i,x) + ^) + v x (B(t,x) + ^)) • V„f = 0, /,U = /o, (1-7) 

for a small parameter e would be a relevant model equation. In equation Q1.7|) . f £ {t, x, v) is at time 
t, the particle density of the plasma in position x and with velocity v. The vector M. € § 2 gives the 
direction of the strong magnetic field. This direction may depend on t and x. the vector AT G § 2 
with M _L Ai is the direction of a strong electric field (VV may also be 0). The vector fields E and 
B are electric and magnetic fields that may contain a self consistent part. We shall call the scaling 
leading to this equation the "Guiding Centre Regime" . 

If now we are interested in understanding what happens close to the plasma boundary in a 
tokamak, we would prefer to use the so called "Finite Larmor Radius Regime" . In other words, we 
would consider an equation of the type: 

^ + V|| • V x f s + ^ • V x f + (E(t,x) + v x ^) • V„/ £ = 0, /,U = fo, (1.8) 

where the meaning of f £ {t, x, v) is the same as above and where V| = (v • .M).M, vj_ = v — vy . 
Lastly, for isotope resonant separation we would be interested in regarding 

| + v ■ VJ e + (E(^,x) + v x (B(i, ^,x) + ^)) • V„f = 0, /|U = /o, (1-9) 

where the electric and magnetic fields E(i, ^^,x) and B(t, ^^,x) oscillate with the same frequency 
l/(2ne) as the cyclotron frequency of the plasma particles. 

We refer to Frenod and Sonnendriicker [TJJ1 HH [HI [14] , Frenod, Raviart and Sonnendriicker 
[9] and Frenod and Watbled [15] where those models are explained and asymptotically analysed. 
In particular, they give indication on how to get the self consistent part of the electric field in 
some examples and under suitable assumptions. Complementary works concerning the asymptotic 
behaviour of those kind of equations are also led in Golse and Saint Raymond [UHH], Saint Raymond 
[32] . Brenier g], Grenier pi] [20], Jabin [21], Schochet [35], Joly, Metivier and Rauch [24]. We also 
refer to mathematical or physical works where similar methods are used: [22l [23l [36] [7] (23 O 

It is relatively easy to see that equations (|1.7[) , (|1.8p and (|1.9[) enter the framework of a singularly 
perturbed convection equation reading: 

du 1 

-£-+a-Vu E + -b- Vu e = 0, u Elf=s =u , (1.10) 
at e 1 

for x £ M. d and t > s. At least formally, the solution u s (t, x) is given by 

u e (t,x)=u (X e (a;x,t)). (1.11) 
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Numerical methods to solve equation (|1.7p . (1 1 - 8[) or (|1.9|) , or similar ones, will be applications 
of generic methods solving (jl.lOp . Since u £ contains high frequency oscillations, using in state any 
well known numerical method to solve (jl.lOp would compel to use a very small time step. 

Yet, in order to relax this constraint, we may follow, at least, two strategies. The first one 
could be based on the work presented in Frenod, Raviart and Sonnendriicker [9] where we proved 
that u £ writes u £ (t, x) = J2i>o £l U l (t, ^T 2 ' x )> an d where we determined the equations satisfied by 
every terms U % . Since those equations are independent of e, the computation of the first terms U % 
could be led using a standard numerical method. Then an approximation of u £ would be given by: 

tt e (t,x) ~ Ei=0 £if/i (*' ^T'X)' for SOme k e N ' 

The second strategy consists in using an approximation X E (t;x, s) ~ J2i=o £ l X 4 (i, x, s), 
of p.5| . for some k £ N, after computing X°,X , . . . ,X fe . Then the expression (jl.lip yields an 
approximation of u e : u £ (t, x) ~ u o(J2i=o £ *X l (s, — ^7^5 x > ^)) ■ This motivates the present work. 

We give now some references concerning perturbed ordinary differential equations. First, we 
think to the Lindstedt-Poincare method explained in Poincare [31] where the steady state periodic 
solutions to a perturbed second order ordinary differential equation is studied. Then the Krylov- 
Bogoliubov-Mitropolsky method, see [3] and [25] , allows to describe the transitory behaviour of the 
solution to a perturbed ordinary differential equation to a periodic solution. 

We also cite the works of Verhulst [H] and Sanders and Verhulst [M] (see also Mickens 
where the Method of Averaging is developed to treat perturbed ordinary differential equations and 
adiabatic invariants in Hamiltonian systems. 

Lastly, we mention the work initiated by Tikhonov [3SJ [391 SO] an( i developed by Vasilieva and 
others (see the review paper of Vasilieva [JT] and the references in it). 

This work is also related to homogenisation methods based on weak—* convergence and two scale 
convergence. The most important ideas about those tools may be found in Tartar [37] , Bensoussan, 
Lions and Papanicolaou [2], Sanchez-Palencia [33] , N'Guetseng [30], Allaire [I] and Frenod, Raviart 
and Sonnendriicker [9]. 

Theorems 

We first justify the expansion |L5j) until order 0. 
Theorem 1.1 Setting 

a°(t,y°) = 77- a°(t,e,y°)d8, (1.12) 



2^ Jo 



with 



a (M,y°) = {V z Z(M;y )} '{afr 9, Z(t, 9; y )) - |^(t, 0; y )} , (1.13) 

under assumptions fl.ty) . fl.3)) with m = and &l-4\ ), for any x £ M. d , s £ R, T £ R and any e > 0, 
the solution X £ (-;x, s) of fZH]) exists on [s, s + T], is unique and the sequence (X £ (-;x, s)) satisfies: 

lim sup |X e (t;x,s) - X°(t, ^^;x,s)| = 0, (1.14) 

£ ^°t£[s : s+T] £ 

| • | standing for the Euclidean norm on W 1 , where X° satisfies 

X°{t,9;K,s) = Z(t,9,Y°(t;x,s)), (1.15) 

and where Y° is the solution to 

rTY° 

^L = a°(t,Y°), Y°( S ;x, S )=x. (1.16) 
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In the Theorem above, V Z Z(£, 9; z) stands for the Jacobian matrix of z i— > Z(t, 9; z). 

In order to set this Theorem and to justify the expansion (| 1 - 5[) for higher orders, using a Van 
der Pol transformation, we define Y £ being such that 

X e (t;x )S ) = Z(f,— ;Y 8 (t,x,*)) = [Z(Y £ )] £ , (1.17) 

e 

where for any function / we write [/(Y £ , . . . , Y*)] £ for f(t, ^ =s -, Y £ (i;x, s), . . . , Y^(t;x, s)) and 
[/(Y°, . . . , Y fc )] e for /(t, ^,Y°(t;x,s),.. . , Y fc (t; x, s)). It is a/easy game to show that 

dY t — s 

^L_£ =5 °(*,l_l Y £ ) = [S (Y £ )] E ; Y £ ( S ;x, S )=x. (1.18) 
Indeed, derivating 17[) and using (|1 . 18[) we get 

= (Ye) ^ + (Ye) ^ + {V,Z(Y £ )} £ {dO(Y e )} £ = [a(X £ )] £ + ib(X £ ), (1.19) 

and we have the initial condition X £ (s; x, s) = x. 

We also use the following notation. For a vector field Z, The i— th component of {V^Zjlx , . . . , x fc }, 
for i — 1, . . . , d, is given by: 

({^ZHx ,* 1 ,... ,**})<= £ J* 2 ' x°...xt (1.20) 

In order to simplify we shall sometimes denote {V^Z}{x°, x°, . . . ,x } by {y^"Z}{x } fe . 
Now for k > we recursively define 

A k (t, 9, y°, . . . , y k ) = ~ | 5 fc (i, a, y°, . . . , y fc ) da - k k {t, y°, . . . , y fc ), (1.21) 

where a° and a are given by ()1.13j) and (| 1 . 12[) and where for k > 1 we have: 

a k (t, 6,y ,..., y k ) = { V„ofi°(i, 0; y°)}{y fe + 0A fc-1 (i, 0, y°, . . . , y fc - x )}+ 
1 fc-i 

-{ V 2 a°(i, 0; y°)} ( £{y^ + 0A J_1 (i, 9, y°, . . . , y^ 1 ) , y k ^ + 0A fc_3 ' -1 (i, 0, y°, . . . , y^ 1 )} 



fc-i 

2 3=1 

1 -{V^<5 o (t,0; y o )}{ y 1 +0A o (t,0, y °)}' £ 



k\ l v 

6(J2{VyiA k -Ht,ey,...y- l )}{&(ty,...,yi)} + ^- — (t,e,y°,...,y*- 1 ) 

3=0 

(1.22) 

1 /"^ 

a fc (t,y°,...,y fc ) = — j o 5 fe (t, 0, y°, . . . , y fe ) d9. (1.23) 



and -, /.27T 



With those notations, we can state the asymptotic expansion of Y £ , 

Y £ = Y° + e(Y x + [0A°] e ) + e 2 (Y 2 + [0A 1 ],) + . . . . (1.24) 
In other words, defining for k > 1, 

Y^ fe = 4r(Y £ - Y° - eCY 1 + [6A°U ^-iry^ 1 + WA^U) - [9A k -%, 

f (1.25) 
= ^(Y*- 1 - Y^ 1 ) - [9A k ~%, 

we have the following Theorem. 
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Theorem 1.2 Under assumptions U.ty) , U.3\) and \1.4\ ) for any x e R d , s G M and T e K, the 
sequences (Y^(-; x, s)), for k = 0, . . . , m, are bounded in L°°([s, s + T]) and we have 

lim sup |Y e (i;x, s) — Y°(i;x, s)| = 0, and lim sup I Y k (t; x, s) - Y k (t; x, s)| = 0, 
te[s,s+T] e ^°te[s,s+T] 



/or fc > 1, where Y° is solution to (1,16\) and where Y k is t/ie solution to 

dY k 



(1-26) 

where Y fc is i/ie solution to 
a fe (t,Y°,...,Y fe ), Y fc (s;x,s) = 0. (1.27) 



As a consequence of this Theorem and of (I1.17p . defining (X^(t; x, s)) by: 

X* = ^(X £ - [X°] £ e^X^ 1 ],) = ^(Xj- 1 - [X^ 1 ],), (1.28) 



for k > 1, we have the rigorous justification of asymptotic expansion (|1.5 



Theorem 1.3 Under assumptions U.S\) . U.3\) and /or any x e M d ; s£l and Tel, £/ie 

sequences (X^'(-; x, s)), /or k = 1, . . . , m, are bounded in L°°([s, s + T}) and we have: 

t — s 

lim sup |X~(t;x,s) - X k (t, ;x,s)|=0, fl.29) 

£ ^°te[s, s +T] e ' 

where X fe writes 

X k (t, 9; x, s) = { V 2 Z(<, 9; Y°(t; x, s))} { Y fc (t; x, s) + flA*" 1 ^, 0, Y°(t; x, s), . . . , Y k ~ 1 (t; x, s))} + 
1 fc-i 

- { V 2 Z(i, 9; Y°(i; x, s))} ( £ { Y J (t; x, s) + 9A^ (t, 9, Y°(t; x, s), . . . , Y J_1 (t; x, s)) , 



Y k -\t; x, s) + ^A^'- 1 (t, 6, Y°(t; x, «),... , Y^ 1 (t; x, s))} 
^{V^Z^, 0; Y°(t; x, s))}{Y 1 (t; x, s) + 0A°(t, 6, Y°(i; x, s))} fc , (1.30) 



wit/i Y fe solution to \1.21^ . 



We end this result list by giving explicit expressions. Using (| 1 . 13[) in (|1.22p we have the following 
expression of a . 

a 1 (t,0,y o ,y 1 ) = {V z Z(i,0;y o )}- 1 {{V K a(t,0,Z(t,0;y o ))}{V,Z( t ,0;y o )}{y 1 +0A o (^0;y )}- 
{^^^y^lly^eA^^^y^l-IV^Z^^y )}^ ^^ )^^^ ^^^ )}}- 

0{{\7 y «A Q (t, 9, y°)}{a°(t, y )} + — (t, 9, y )) . 

(1-31) 

On another hand, if b is linear and independent of t, a 2 is given by 



a 2 



(t 7 9y,y\y 2 ) = {X7 z Z(9;y°)} 1 { { V,a(i, 9, Z(9; y ))} { \7 Z Z(9; y )} {y 2 + 9A 1 (t, 9; y°, y 1 )} 
i{V 2 2 a(t,0,Z(0;y o ))}{{V 2 Z(0;y°)}{y 1 +0A o (t,0;y )}} 2 - 

0({V,o A^t, 0, y°, y 1 )}^, y )} + {V„i A^t, 0, y°, y 1 )}^, y°, y 1 )} + -^(t, 9, y°, y 



(1.32) 



The paper is organised as follow: in the second section we briefly prove the Theorems. Section[3]is 
then devoted to applications to models describing magnetic confinement fusion and isotope resonant 
separation experiments. Among those applications, one involves a non uniform strong magnetic 
field, in this case the computations are led using Maple. 

Remark 1.1 With very little changes we can apply the previous results to the case when b = 
b(i, 9, x) also depends on 9, as soon as the regularity of b is enough and the assumption (|1.4|) is 
realized. 



2 Proof of the Theorems 

2.1 Sketch of the proof of Theorem 11.21 

The key point to prove the results of this paper is Theorem 11.21 in the same spirit as in Schochet 
[55] , This result is classical and known as Single Phase Averaging. We do not give a detailed proof 
of it but only a formal sketch of it using expansion (| 1 . 24[) of Y £ . 

First, because of the definition (|1.13p of 5° and the assumptions (|1.2|) - (jl.4p we deduce that the 
function (i, y°) i— > a°(t, ^-,y ) is regular enough to ensure that, for any e > 0, Y e (-;x,s) exists, is 
unique on [s, s + T] and remains in a bounded set of M. d independent of e. 
Then expanding &°{Y e ), using Y s = Y° + J2j>i ^( YJ + [Q&'^e), we obtain 

&°(Y E ) = &°(Y°) + e^oc^Y )}^ 1 + {9A%} + e 2 ({V,oa°(Y°)}{Y 2 + [9A%} + 
1 



{VJUWHY 1 + [6k\Yj +---+e k (jV,oa°(Y°)}{Y fc + [flA^ 1 ],}- 
1 fc-i 

-{V y 2 o5 (Y°)}(^{Y^ + [0A-''- 1 ],, Y fe ^ + [9A k -J- 1 } e }) + 



±{V y k oa°(Y°)}{Y' + l9AX} k 



(2.1) 

Plugging the expansion (|1.24[) in the dynamical system (| 1 . 1 8[) and identifying the terms of the same 
order, we have 

9Y° 3[0A°]_ 

~dT + ^9-~ a {t ' e ' Y } ' (2 - 2) 

at the order —1, and at the order k — 1 for k > 1 we have 



8Y k d[9A k ] _ . A 



dt 09 



a K (t,9,Y\...Y K ). (2.3) 



Then equation (fTTTCj> for Y°, definition (TOT]) of A fc and equation (Tj~2T|) for Y fc follow, ending the 
sketch of the proof of Theorem 11.21 ■ 

2.2 Proof of Theorem fTTTI 

Existence and uniqueness of X e are consequences of existence and uniqueness of Y e . Then the 
regularity of Z and Theorem 11.21 allows to pass to the limit in (| 1 . 1 7[) in order to deduce (|1.14[) and 
(|1.16p and thus prove the Theorem. ■ 
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2.3 Proof of Theorem [TT31 

From Theorem II. 2 [ we deduce 

Y E = Y° + £ (Y X + [0A°] e ) + • • • + e*- 1 ^* -1 + [AA^ 2 ],) + £ fe (Y* + [9A k - 1 ] e ), (2.4) 
and then, from (| 1 . 1 T[) we deduce 



X £ = [Z(Y e )] e = [Z(Y°)] e +e{V,Z(Y £ )} £ {Y 1 + {6A°} £ } + e 2 (jV z Z(Y £ )} e {Y 2 + [0A 1 ] e } + 

i{V 2 Z(Y £ )} £ {Y 1 + [9A\} 2 ^j +■■■ + £ fc ({V,Z(Y e )} £ {Y* + [0A fc - 1 ] £ } + 
1 fc-i 

-{V 2 Z(Y £ )} £ (]T{Y^ + [^•- 1 ] E) Y fc -'- + [9A k -^] e }) + 

■■■ + ^{V^Y )}^ 1 + [0A°] e } k> j + o(e k ). (2.5) 

Comparing this expansion with 

X £ = [X ] e + e[X 1 ] e + .» + e fc - 1 [X fc - 1 ] e +e fc Xj J (2.6) 
obtained from (|1.28[) and making the process e — > give finally Theorem II. 31 ■ 

3 Application to the gyrokinetic plasma 

The dynamical systems associated with equations 111. 7)1 . OH} and (|1.9| are in the form of (|1.5p . i.e.: 



* \V t/ ) e e \V,(s;x,v,.,) 

with variable (x, v) G R 3 x R 3 in place of x and with ad-hoc fields a and b. Then we can apply our 
result saying that the solution (X e (£; x, v, s), V £ (t; x, v, s)) can be expanded as 

XJt; x, v, s) = X°(i, - — -; x, v, s) + eXVi, - — -; x, v, s) + e 2 X 2 (i, - — -; x, v, s) + . . . , 

t-s t-s t-s ( 3 - 2 ) 
VJt; x, v, s) = V°(t, -;x, v, s) + eV 1 ^, ; x, v, s) + e 2 V 2 (i, ; x, v, s) + . . . , 

e e e 

and that this expansion may be justified until any order if the regularity of the fields is enough. 

In this section, we first lead the computations in the case of Isotope Resonant Separation Regime 
with the restriction that M. is constant and B(t, 0, x) = until the order 1. Then, with the same 
restrictions, we treat the case of the Guiding Centre Regime until the order 2 and we give the result 
for the Finite Larmor Radius Regime until the order 0. The forth example concerns the Guiding 
Centre Regime with a variable strong magnetic field and a constant B until the order 1. 

From the physical point of view, the last example is relevant to understand the behaviour of a 
plasma in a tokamak. This example shows that the generic computations made before, coupled with 
the use of Maple, enable to deduce the result relatively comfortably. 
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3.1 Isotope Resonant Separation Regime with constant strong magnetic 
field 

In the case of Isotope Resonance Separation Regime, i.e. of equation p.9[) . a and b are: 



a(t,0,x,v) 



E(i,0,x) + v x B(t,0,x) 



, b(t,x,v)=b(v) 





v x M 



(3.3) 



For simplicity, we restrict to the case when B(i, 0, x) = and when M. = e\ is a constant vector, 
(ei,e2,e3) being the frame of M 3 . Then Z(i, 0; z, w) and {V 2jlu Z(i, 0; z, w)} _1 are given by 



Z(i,0;z,w) 



, {V ZiU ,Z(i,0;z,w)} 1 = 



z 

#(0)w y 

where i?(0) is the matrix of the rotation of angle —0 around Ai, 



I 
R{-9) 



(3.4) 










\ 


i?(0) = 





cos 


sin0 




1° 


— sin 


cos 9 j 



(3.5) 



We have the following result. 

Theorem 3.1 If~E(t, 0,x) is C,J(lxRx M. 3 ) and 2n— periodic in 9, the first term of the expansion 
of the solution (X e (i; x, v, s), V E (i;x, v, s)) £o 

^=V £ , ^=E(t,— ,X £ ) + iv e xM, X £ ( S ;x,v, S ) =x, V e (s; x, v, s) = v, (3.6) 

at at e e 

is given by 

X°(t,0;x,v,s) = Y°(i;x,v,s), V°(t, 0; x, v, s) = i?(0)U°(i; x, v, s), (3.7) 
where (Y°(t; x, v, s), U°(i; x, v, ,s)) is solution to 

V = U "' ~dT = 2^l ^(-W^» Y °)^> Y°( S ;x,v, S )=x, U (s;x,v, S )=v. (3.8) 
Proof. In view of (|1.13[) and (|1.12p . here we have 



a°M,y ,u°) = 



i?(0)u° 
«(-0)E(t,0,y°) / 



, a°(t,y°,u ) 



2tt 



(3.9) 



Then the proof of the Theorem is straightforward. 



In order to obtain the system satisfied by the second term ( Y 1 , U 1 ) of the expansion we notice that 
U(9) = -i?(f+0) + i?(§) is such that R{a) da = 9P + K(9), with P the matrix of the orthogonal 
projection onto A4. We have 



(3.10) 





^0 









(l 





o\ 


n(B) = 


sin 6 


1 — cos 


9 


and P = 













^0 cos - 


-1 sin0 


J 




[° 








Then, from (|1.21|) . we get: 

<?A°M,y ,u°) 



TZ{9)u° 

( l^-^Jdl)(R(-a)E(t,a,y )) 



(3.11) 



where we denote f® f(a) da - ^ f{a) da by ( J® da - ^ ^ da) (f(a)). We also get 



-1U o i u W)(u X +( [da-^-[da)(R(-a)E(t,a,y 
a 1 (t,9,y ,u?,y 1 ,u 1 )= \ Wo 27r J Q /V 

\ R(-0)V x E(t, 9, y°)(y 1 + ft(0)u°) 

W ° 2Wo i?(-o-) V x E(i, a, y°)ujj 







(t,«7,y°) 



(3.12) 



Then using ( / dcr - ^ J Q da) (i?(<r)) = K{0) and i /„ K{9) d9 = i?(f ) - P, we obtain 
..o ..o ..1 „u _ 

\ 



u S + - y ( y^ ^ - y*,) - . )E( , „, ^ * 
i /" z7r i /" 27r 

VW fi (-^) v - E (*.^y )^)y +(^y R(-0W x E(t,6,y o )1l(9)d9)u 
' (R(~)-P)(±j\-<;Mt,<;,y )d<;) ^ 



2tt , /.e 



\ V 27r 

/ 

1 

Hence we can state the following Theorem. 



1 r / r o r \ / \ 



2tt 



2tt, 



<9E 



do- - — / da ){R(-a) — (t , a, y° ) otf 



at 



(3.13) 
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Theorem 3.2 Z/E(i,0,x) isC 2 (JR xRxl ) and 2tt — periodic in 9, the second term of the expansion 
of the solution (X E (t;x, v, s), V e (t;x, v, s)) to 6')) is given by 

X 1 (i, 0; x, v, s) = Y 1 (t; x, v, s) + K(9)XJ° (t; x, v, s), 

V x (t,0;x,v,«) = i?(6»)U 1 (t;x,v,s) +i?(<9)( da - — da)(R(-a)~E(t,a,Y°(t;x,v,s))), 

Jo 2ttj 



(3.14) 
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where (Y^U 1 ) is solution to 
dY 



dt 



/ R(~9)V x E{t,9,Y°)d9)Y 1 + (— R(-9)W x E(t, 9, Y°)1l(9)d9)U a - 

Jo . 27T J Q 



dt 2n J Q 

( tJ" (ff ~ UP) (fl(-)V,E(^,Y»))<i e )U5- 
Y 1 {s;x,v,s) = 0, U x (s;x,v,s) =0. 

(3.15) 

3.2 Guiding Centre Regime with constant strong magnetic field 

In the case of the Guiding Centre Regime, i.e. of equation (jTTTJ) , we have 

/ v \ / 

a(t, 9, x, v) = a(t, x, v) = , and b(t, x, v) = b(v) = 

\E(i,x) +v x B(t,xW \A/" + vx.M 

(3.16) 

As previously, we make the restriction B(f, 0, x) = W = and .M = ei. Since this situation 
is similar to the previous one with the only difference that E(t, x) does not depend on 9 we can 
directly deduce the following Theorem. 

Theorem 3.3 7/E(t,x) is C^(R x M 3 ), the first and second terms of the expansion \3.2\) of the 
solution (X E (t; x, v, s), V e (t;x, v, s)) io 

^=V e , ^=E(a £ ) + -V £ x^, X e ( S ;x,v,s) =x, V e (s; x, v, s) = v, (3.17) 
at at e 

are given by 

X°(t, 0; x, v, s) = Y°(t; x, v, s) , V°(t, 0; x, v, s) = P(0)U°(i; x, v, s), (3.18) 

and 

X 1 (t, 0; x, v, s) = Y x (i; x, v, s) + ft(0)U°(t; x, v, s), 
V 1 (i, 0; x, v, s) = R(9)JJ 1 (t; x, v, s) + ft(0)E(t, Y° (t; x, v, «)), 
where (Y°(i; x, v, s), U°(t; x, v, s)) is solution to 
dY° d\l° 

US, — =E|,(t,Y°), Y°(s;x,v, S )=x, U°(s; x, v, s) = v, (3.20) 



(3.19) 



dt 11 dt 

and where (Y^U 1 ) is solution to 



dY 1 



Un + (R{ J) - P)E(t, Y°) = U^ + E(<, Y°) x M, 



dt 11 v '2 

^- = P VrE(t, Y°)Y 1 + itr((J - P)V*E(i, Y°))(P(-|) - P)U° + 

^ r ((^(-|) - ^)V,E(t, Y )) (/ - P)U° - (P(-|) - P)V x E(i, Y°)U° - ( 3 - 21 ) 

W-f)-P)f(*,V0), 

Y x (s;x,v,s) = 0, U^s^v.s) = 0. 
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In the computations leading to this Theorem, we use, among other formula, 

(Jda-—Jda)(R(-a))=-K(-9), —J -K(-9) d9 = R(--) - P, -R(9)K(-9) = K(9), 

(3.22) 

and 



— J R(-e)V x E(t,Y°)n(0)d6 = -tr((I - P)V x B(t,Y°))(R(-^) - P) + 

itr((P(-|)-P)V x E(t,Y°))(I-P). 

Now we turn to the characterisation of (X 2 , V 2 ). For this purpose we notice that here 
a 1 (t,e,y ,u°,y 1 ,u 1 ) = 



(3.23) 



P^u 1 + ft(0)E(i,y°) 



<9E 



R(-9)V x E(t, yV + R(-9)V x E(t, y°)K(9)u + K(-9)V x E(t, y°)u° + H(-0) — (t, y°) , 

11 at 



In order to get now 9 A 1 we need first to compute 

( jda - — J da) (K(a)) =1- R(9) , ( J da - — J da) (K(-a) ) = I - R(-9). 



(3.24) 
(3.25) 



Secondly, 



Jda - ^Jda) (P(-<7)V x E(i, y a )U(a)) 



PV a E(i, y°)(J - R(9)) + (I - P)V,E(i, y° 



A) 



v° 





sin 2 9 



\ 

sin 26* 



sin ( 



sin 26» 



sin 2 (9 



(p(-|)-P)V x E(t,y° 



sin 





sin26> 

+ cos 6» - 1 



1 — cos 9 — 



sin 



- sin 26» 



which also reads 

(.27T 



(3.26) 



f f) f 

( da-— da)(R(-a)V x E(t,y°)K(a)) = PV x V(t,y°)(I - R(6)) 
Jo Z7T Jo 



\{ - n{-9) + P(|) - p)v x E(t, y °)(^) + p(|) - p) + §(-«(— - P)v x nt,y°)(R(^) - p). 

(3.27) 

Hence integrating (|3.24p we have 



9A 1 (t,9,y°,u° 1 y\u 1 ) = 



K^u 1 + (J - R(0))E(t,y°) 



-tt(-0)V x E(t, y^y 1 - (I - R(-9)) (v a E(t, y°)u{| + ^(t, y°) 
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PV x E(t, y°)(I - JZ(ff)) + - ( - + R{\) - P)V X E(«, y°) (11(6) + R(-) ~ P)j u° y 



o 

7T I • (3-28) 

P)V x E(t,y°)(i*(|)-P)u / 



In order to obtain the expression of a , we need to compute 
{V, lU ,Z(t, 0; y^u )}"^^,^. > ; y°. u °))}{ V,,„Z(t, 9; y°, u°)} 



P(6») 
fl(-0)V x E(i,y°) / V\u 



0A 1 (f,0,y°,u ( \y 1 ,u 1 ) 



/ i?(0)u 2 +^)V K E(i,y°)y 1 - - /) (v a E(t, y°)u|j + ^(*,y°)) + 

(pV x E(i,y°)(J - P(0)) + - P)V*E(t,y°)(ft(0) + iJ(f ) - P) + 

|(£(0-f)-P)V x E(i,y o )(ii(f)-P))u o 



V 



fl(-0)V z E(t, y°) (y 2 + ft^u 1 + (I - H(0))E(t, y°) 
{V 2 a(i, 0, Z(t, 0; y°, u ))} J {v z ,„Z( f , 0; y°, u°)}{( y J + 0A°(i, 0, y°, u )} 



(3.29) 



(3.30) 



v { V 2 E(t, y )} 2 ({y 1 , y 1 } + 2{y\ ^)u } + {TZ(9)u°, Tl(9)u°} 



A direct but heavy computation also leads the derivatives V y o iU o A 1 (t, 9, y , u , y 1 , u 1 ) a°(t, y , u ), 

V y i, u iA (t, 9 7 y°, u°, y 1 , u 1 ) a 1 ^, y°, u°, y\ u 1 ), and ^-(t, 0, y°, u°, y 1 , u 1 ). 

Then in view of (11.321) . we get a 2 , which, integrating with respect to 9 leads to a 2 and to the 
equation satisfied by (X 2 , V 2 ). Then we have the following Theorem. 

Theorem 3.4 If E(t, x) is Cf(M x R 3 ), the third term of the expansion \3. S\) of the solution 
(X e (i;x, v, s), V E (t;x, v, s)) to {3.17^ is given by 

X 2 {t, 9; x, v, s) = Y 2 (i; x, v, s) + ^(^)U 1 + (/ - P(0))E(t, Y°(t; x, v, «)), 
V 2 (i, 0; x, v, s) = i?((9)U 2 (t; x, v, s) + ft(0)V x E(t, Y°(t; x, v, s))Y 1 (i; x, v, s) - 
(R(9) - J)V B E(t, Y°(t; x, v, s) (u«(t; x, v, s) + — (f, Y°(t; x, v, a))) + 

(pV x E(t, Y°(i;x,v, S ))(/ - + i(i?(|) - P)V x E(t, Y°(i;x,v, S ))(ft(0) + - P) + 

l -(R(9 - |) - P) V a E(t, Y°(t; x, v, s))(P(^) - P))u°(t; x, v, s). 

(3.31) 
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Moreover (Y 2 ,U 2 ) is solution to 



Uf + - P)V x E(t, Y^Y 1 + (7- P)(v x E(i, Y°)U|j + ^(t, Y°) 



(PV x E(t, Y°)(7 - P) + (P(|) - P)V x E(i, Y°)(P(|) - P))U°- (3.32) 
itr((7 - P)V x E(i, Y )) (7 - P))U° - ~tr((P(-|) - P)V x E(t, Y°))(P(|) - P)U°, 

^ = PV x E(i,Y°)Y 2 + 

~tr((7 - P)V x E(i, Y )) (P(-|) - PjU 1 + itr((P(-J) - P)V x E(i, Y°))(7 - P)IP + 
PV x E(t, Y°)E(t, Y°) - (PV x E(t, Y°)P + ^tr((7 - P)V x E(t, Y )) (7 - P)+ 
^tr((P(-|) - P)V x E(i, Y°))(P(-|) - P))e(*, Y°)+ 



U V 2 E(t, Y )} 2 ({Y 1 , Y 1 } + 2{Y 1 , U|}} + {(7 - P)U°, (7 - P)U }- 



{(P(-f)-P)U°,(P(-|)-P)U } 

(^- P(-f )){V|E(t, Y°)}{Y 1 ,U|}} + (7 - P)({V 2 E(i, Y°)}{UU,U|f} + ^(t, Y°)Ujf)- 
P{V 2 E(t,Y°)}{(7-P)U° ! U«}- 

i(7-P){V 2 E(t,Y°)}{(7-P)U°,U{J}- ^(^(-^-^{V^Y^H^^-^U ,^}- 
- (7 - P) (V x E(t, Y )) + PV x E(t, Y°)(7 - P)+ 



~(7 - P)V x E(t, Y°)(7 - P) + -(«(--) - P)V x E(t, Y°)(i?(|) - P))jZ\\(t, Y°)+ 
(P - P(-|))V x E(t, Y°)U} + (P - P(-|))V E E(t,y°)(P(|) - P)E(i, Y°) + 

(p -^f )) ^ (^ ' YO)Yl - (/ - p K^ (^ ' YO)u ^^ (^ ' YO 0- p ^^ (t ' YO)(/ - p) - 

i (/ - p) ^lr (< ' Y ° )(/ - p) + zW-$ p )^r^ Y °^ R fy p )) u °- ( 3 - 33 ) 

3.3 Finite Larmor Radius Regime with constant strong magnetic field 

In the case of Finite Larmor Radius Regime l|1.8[) with M. = ei , we have 

a(i,0,x,v) = I V " and b(<,x,v) = b(v) = [ V± I, (3.34) 

(z + K(9)w\ , (i K(-9)\ 

Z t,M,w)= , {V, lto Z(i,0;x,v)} = M. (3.35) 

y R(0)vr ) \0 R( 

with R and ft defined by ([33]) and (|3~TTJ)) . 

Here we give the result for this case only for the order 0. We have 

a°(i,0,y°,u°) = 11 K ' K / n h (3-36) 

P(-6»)E(y° + ft(0)u o ,t) 
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and 



a°(t,y ,u°) = 

V ^J 

Hence we have the following Theorem. 



ft(-0)E(y° +H(0)u°,t) d9 
R(-9)E(y° + TZ{9)u°,t) d9 



\ 



(3.37) 



Theorem 3.5 If we assume that E(f, x) is C^(R x R 3 ), the first term of the expansion 113. 2\) of the 
solution (X E (t; x, v, s), V E (t; x, v, s)) to 

^=V e ||+iv e± , ^ = E(t,X E ) + -V £ xM, X e ( S ;x,v,s)=x, V e ( S ; x, v, s) = v, (3.38) 
at e etc e 

is given by 

X°(i,0;x,v,s) = Y (t;x,v,s)+K(«)U (t;x,v,s), V° (i, 0; x, v, s) = R(6)XJ°(t; x, v, s), (3.39) 
where (Y (t; x, v, s), U°(t; x, v, s)) is solution to 



dY° 

dt 



1 f 2Tr 

Uy + 7^ / ft(-0)E(Y° + 7£(0)U°, t) d0 , 

11 27r Jo 



o!U° 1 



dt 2tt 



/■27T 

/ i?(-0)E( Y° + TZ(9)V° , t) d0, 
Jo 

(3.40) 



with the initial conditions 



Y°(s; x, v, s) = x , U°(s; x, v, s) = v. 



(3.41) 



3.4 Guiding Centre Regime with variable strong magnetic field 



Here we study a situation representative of what happens inside a tokamak, i.e., the Guiding Centre 
Regime with a variable M, with J\f = and B = e3. In other words we consider the following 
system: 

' /X V,, ^ = E(i,X E ) + V £ xe 3 + -V f xW(X E ), X e (s; x, v, s) = x , V e (s; x, v, s) = v, 



dt 
where 



dt 



M(x) 



\lx\ z + x 2 2 



(-x 2 \ 

Xl 

V J 



p T (x)ei, with p(x) = 



\/xi 2 + x 2 2 



where x = (xi,X2,xs) in the frame (ei, 62,63) of R 3 . In this case 



a(t, 9, x, v) = a(t, x, v) 

' E(i,x) + v x e 3/ 
and thus, R{9) being defined by (|3.5p . 



, and b(i, x, v) = b(x, v) 



-x% x\ 0^ 

-X\ X 2 

1 





v x M(x) 



(3.42) 
(3.43) 



(3.44) 



Z(t,0;z,w) 



T (z)fl(0)p(zW 



, where p T Rp 



( Z2 2 +Zi 2 cos(#) ^1^2 (cos(#) — 1) z\ sin(#) \ 

2l 2 + 22 2 Xi 2 +Z 2 2 VZ1 2 +Z 2 2 



Zl 2 + Z2 2 

jSl sin(0) 

Vzi 2 +Z2 2 



Zi 2 +Z 2 2 %/zi 2 +Z 2 2 

z 2 sin(fl) 



Vzi 2 +Z2 2 



and we have the following Theorem. 



cos(0) J 

(3.45) 
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Theorem 3.6 J/E(t,x) is C 2 (Rx]R 3 ) 7 the first term of the expansion VS.ty of the solution (X e (i;x,v, s), 
V £ (i;x, v, s)) to Jff.^£| ) is gwen 6y 

X°(t, 0; x, v, s) = Y° (*; x, v, s) , V° (t, 9; x, v, s) = p T (x)i?(%(x)U° (t; x, v, s), (3.46) 

where (Y°(i; x, v, s), U (t; x, v, s)) is solution to 

d\J° 



dY° 

^ = A(Y°)U°, 



wii/i .A(y) = 



= /3(Y U , U u ) + A(Y )Eii (t, Y") + U u x A(Y u )e 3 



1 



/ y 2 2 



yi 2 + y2 z 



-2/12/2 



\ 



V 



-2/12/2 yi 2 o 
o o oy 



, and /3(y,u) 



( {V2U\ - yiU2)u 2 \ 
Vi 2 + V2 2 

(VlU2 ~ V2Ul)U 1 



\ 



yi 2 + y2 2 
o 



(3.47) 



/ 



Y°(s;x, v, s) = x, , U°(s; x, v, s) = v. 
TTie ierm X 1 (t, 0; x, v, s) is given by 

X\(t,0;x,v,a) = i(Y° 1 fi(a*(0) - 1)U° 3 + Y° 1 (Y° 1 U° 1 + Y° 2 U° 2 ) sin(0) + Y 1 ^ 2 ) , 
X 1 ^*, ^; x, v, s) = (y° 2 Q(cos(0) - 1)U° 3 + Y° 2 (Y° 1 U° 1 + Y° 2 U° 2 ) sin(0) + Y 1 ^ 2 ) , 
XSft 0; x, v, s) = I ^(— Y°iU°x - Y° 2 U° 2 )(cos(0) - 1) + Y\ 3 ft + sin(f?)U° 3 ft) , 

(3.48) 

where ft = \/Y i 2 + Y° 2 2 and where Y 1 (t;x, v, s) is solution to 
dY 1 1 



dt 



ft 2 Y° 2 U° 2 + 2Y 1 Y° 2 (-U 1 Y° 2 + U° 2 Y°i) ) Y 1 ! + (2U° 1 Y° 2 - U^Y ^ 



2Y° 2 '(-U° 1 Y° 2 + U° 2 Y°i)j Y 1 2 + Y° 2 U\ft 2 - Y 1 Y 2 n 2 V 1 2 - n 3 Y\B 3 - Y° 2 U° 3 ft 3 - 

2Y° 2 U 3 (-U 1 Y 2 + U 2 Y° 1 )^ /n\ 



dt 



((2U° 2 Y°i - lAY^f! 2 - 2Y° 1 2 (-U° 1 Y° 2 + U° 2 Y°i)) Y 1 ! + ( - n 2 Y°iU°i- 
2Y 1 Y° 2 (-U\Y° 2 + lAY !)) Y 1 , - Y° 1 Y° 2 f7 2 U 1 1 + Y^U 1 ^ 2 - 

ft 3 Y° 2 E 3 + £l 3 Y°iU° 3 + 2Y Q 1 U 3 (-U° 1 Y 2 + UWiJfi') / ft 4 , 



dY 1 



dt 



= ft 2 Y°iEi + ft 2 Y° 2 E 2 + ( U^Y ! + UY - U°! Y° 2 + U° 2 ) ft 



(yW + y° 2 u 2 ) j /ft 



Y 1 (s;x,v,s)=0, 

(3.49) 

with E being evaluated in Y . The result concerning V 1 (t, 9; x, v, s) is given in the avvendix HI 
This Theorem is the consequence of the following computations. Setting A = p T Rp we have 

{V z , w Z{t,9-z,w)} = 1 ° and {V 2 ,„,Z(<,0;z,w)} _1 = 1 ). (3.50) 

\V,(Aw) A/ ' U T V 2 (Aw) A T / 
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Then, since in this example = 0, and 



Aw 

a(t,Z(M;a,w))= | |, (3.51) 

E + Aw x e 3 



applying (|1.12p and (|1.13|) we get 



a°(f,0,z,w) = I AW I, (3.52) 

- A T V z (Aw) (Aw) + A T E + w x A T e 3 



and since J v Add = J Q W A T dO — 2itA, and L " — A T V z (Aw)(Aw) = /3, we finally get the first part 
of the Theorem. 

The second part of the Theorem is obtained following the computation program described in the 
previous sections using Maple. 

A Appendix : V 1 for the Guiding Centre Regime with vari- 
able strong magnetic field 

The computations of this appendix have been realized using Maple5.5, on computers of Medicis 
Centre, Ecole Poly technique, Palaiseau, France. 

The electric field E and its derivatives being evaluated in Y° , and Q standing for -\/(Y i) 2 + (Y°2) 2 , 
the term V 1 (t, 9; x, v, s) of the expansion (|3.2p of the solution to (13.42|) is given by : 

V 1 1 (i;x,v, S ) = M ((2Y° 1 U° 1 + Y° 2 U° 2 )fi 2 - 2Y° 1 2 (Y° 1 U° 1 + Y° 2 U° 2 )) Y\+ 

(ViUW - 2Y° 1 Y 2 (Y iU° 1 + Y° 2 U° 2 )) Y x 2 + yVuW + Y°iY 2 n 8 U 1 2 + fi 3 Y°iE 3 + 
O 3 Y 2 U°3 + Y° 2 U 3 (-U 1 Y 2 + Y°iU°a)n)oo6(0) + f(-" 3 U° 3 + fiU° 3 Y° 1 2 )Y 1 1 + 



Y iU 3 Y 2 nY 1 2 - Y° 1 n 3 U 1 3 + f7 2 Y°i 2 E! + fi 2 Y°iE 2 Y° 2 + U° 2 ft 4 + 



U 2 (-U°iY° 2 + Y O iU o 2 )fi 2 Jsin(0) + ( - fi 2 Y° 2 U° 2 + 2Y 1 Y° 2 (-U 1 Y° 2 + Y°iU° 2 )j Y 1 i + 
(2U\Y° 2 - Y°iU° 2 )n 2 + 2Y° 2 VU° 1 Y° 2 + Y°iU° 2 )^ Y 1 2 + Y° 2 2 uW - Y'^yWu 1 ,- 

f7 3 Y°iE 3 - n 3 Y° 2 U° 3 - Y 2 U° 3 (-U 1 Y° 2 + Y°iU° 2 )n ] / fi 4 , (A.l) 



V\(t;x,v,s) = (lAY 2 ft 2 - 2Y 1 Y° 2 (Y° 1 U () 1 + Y° 2 U° 2 ))Y 1 1 + (Y'^U ! + 2Y° 2 U° 2 )n 



2Y° 2 2 (Y° 1 U° 1 + Y° 2 U° 2 )J Y 1 , + Y 1 Y° 2 Q 2 \J\ + Y^U 1 ^ 2 + fi 3 Y° 2 E 3 - Q 3 Y X U 3 - 

Y 1 U 3 (~U 1 Y° 2 + Y°iU° 2 )nJcos(6>) + fY iU 3 Y 2 OY\ + (-fi 3 U° 3 + nU° 3 Y° 2 2 )Y 1 2 - 

Y° 2 fi 3 U 1 3 + fi 2 Y° 1 Y° 2 Ei + f7 2 Y° 2 2 E 2 - U°ifi 4 - U°i(-U°iY° 2 + Y°iU o 2 )fi 2 ^sin(0) + 

^2Y 1 U 2 -U 1 Y 2 )fi 2 -2Y 1 2 (~U 1 Y 2 +Y 1 U 2 ))Y 1 1 + 
(-n 2 Y°iU"i - 2Y 1 Y t, 2 (-U°iY 2 + Y°iU° 2 ))Y 1 2 - Y°iY° 2 n 2 U :L i + yVu 1 ^ 2 - 

fi 3 Y° 2 E 3 + ^Y ^ , + Y°iU 3 (-U 1 Y° 2 + Y°iU° 2 )n ) / fi 4 , (A.2) 



1G 



V 1 s(t; x, v, a) = ^^U^ - n 2 Y°iEi - fi 2 Y° 2 E 2 + (U"iY () 2 - Y°iU° 2 )fi 2 - 

(-U !Y° 2 +Y° 1 U° 2 ) 2 )cos(0) + ((lAfi 2 - Y° 1 (Y° 1 U° 1 + Y° 2 U° 2 ))Y 1 1 + 
(U° 2 fi 2 - Y° 2 (Y°iU°i + Y° 2 U° 2 ))Y 1 2 + Y° 2 n 2 U 1 2 + 3 E 3 + Y o ifi 2 U 1 i > )sin(0)+ 



a 2 Y°iEi + n 2 Y° 2 E 2 + (U°i 2 + Y°iU° 2 - U°iY° 2 + U° 2 2 )a 2 - (Y !U°i + Y° 2 U° 2 ) 2 J J fi 3 , (A.3) 
where U 1 (t; x, v, s)) is solution to: 

^ = 1/4 (("4 Y° 2 Q 8 Y°! §5| +4Y° 2 2 O 8 |§i ) Y 1 3 + ( - 4f2 9 Y° 2 + (-yV U° 3 +4Y° 2 U° 2 Y% - 4U\ 

Y° 2 2 ) Q 7 + (-YV U° 1 -2Y° 1 Y° 2 3 U°i+2 Y° 2 4 U° 1 +2Y° 1 2 Y° 2 2 U° 2 +Y° 2 2 Y° i 2 U° 3 -U° 2 Y° 2 Y\ 3 - 2 yV 
U°i Y° 2 2 + YV U°3 Y° 2 - 4 U° 2 Y° 2 3 Y°x) +Y ! (3 yV Y° 2 2 + Y° 2 4 + Y°i 4 + Y° 2 Y°i 3 ) (Y x U°i + Y° 2 
U° 2 ) fi 3 ) E 3 +(4 n 8 Y° 2 U° 2 — Y°i 3 f2 6 U° 3 +Y° 2 Y°i 2 (-Y° 2 2 U° 2 - U°i Y° 2 Y°i + Y\ 2 U° 3 + U° 3 Y°i Y° 2 ) 

O 4 +Y ! 4 Y° 2 (Y^ + YM (Y<\ U° 1+ Y° 2 U° 2 ) O 2 ) U 1 1 + (-f2 7 Y<\ U° 3 +Y° 1 (8 Y°i 2 U° 3 -Y° 2 U° 2 Y° 1 — U° 1 
Y° 2 Y°! + 5 U° 3 Y<\ Y° 2 + 2 Y° 2 2 U° 3 ) n 5 - Y°i 2 (8 Y<\ U° 3 Y° 2 2 - Y<\ Y° 2 2 U°i + 7 Y°i 2 Y° 2 U° 3 - Y<\ Y° 2 2 
U° 2 + 6 Y° 2 3 U° 3 + 7 Y°i 3 U° 3 - YV Y° 2 U°i - Y° 2 3 U° 2 ) O 3 ) Y 1 ! 2 + (yV Y° 2 5 U° 2 3 - 8 uV U° 3 Y° 2 7 + 7 



UV U° 3 Y°i' + Y° 2 2 YV U , 3 + 6 U° 3 YV U° 2 2 + Y°i 6 U°i 3 Y° 2 + Y° 2 4 Y°i 3 U° 2 3 + 5 U° 2 2 Y°i 5 Y° 2 2 
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